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I. INTRODUCTION 

Finite volume/finite size effects play or may play an 
important role in systems ranging from the physics of 
ultracold atom clouds or optical lattices over condensed 
matter systems, multi-layer systems to heavy-ion colli¬ 
sions where the interaction region is relatively sharply 
bounded, though expanding. On the theoretical side, lat¬ 
tice methods or other theoretical approaches with space- 
time coarse graining are subject to finite volume effects 
due to finite number of lattice or grid points. 

In an Euclidean field theory finite volume effects are 
described similarly to that of thermal fluctuations, hence 
being related to a wealth of interesting phenomena and 
showing characteristic scaling laws. In particular, sys¬ 
tems with a second order phase transition lose this prop¬ 
erty in a finite volume and show finite size scaling. Promi¬ 
nent questions in this context concern the approach to 
the infinite volume limit and the continuum limit in lat¬ 
tice theories, the characteristics of the finite volume scal¬ 
ing and the appropriate extraction of thermodynamical 
observables in a given approach. 

It is well-known from systems at finite temperature 
that thermal effects both show a characteristic exponen¬ 
tial decay with the masses of the system at hand, as 
well as a characteristic polynomial decay with momenta 
due to thermal contact terms. The latter contact terms 
also play a pivotal role in thermodynamical relations. 
Moreover, the strength of such contribution is only com¬ 
puted with the correct exponential thermal decay. It is 
this peculiar combination of exponential and polynomial 
suppression that requires a quantitative control in non- 
perturbative approaches. 

In the present work, we discuss finite temperature and 
finite volume effects within the functional renormalisa¬ 
tion group (FRG) approach, see e.g. [TH5]. Here, we ex¬ 
tend the discussions and approach set up in Refs. EHin] 
with a special emphasis on the interplay of the scaling 
with temperature T, volume (box length L), and with 
the infrared cutoff scale k in the FRG approach. This 


is worked out at the example of a real (/)^-theory. Only 
if all these scalings are taken into account quantitatively 
one can expect to have access, in particular, to thermo¬ 
dynamical observables and the equation of state. 

A related functional approach is given by Dyson- 
Schwinger equations (DSEs) at finite temperature and a 
finite volume, see e.g. pTlllJ) . if the infrared cutoff scale 
k is substituted with the ultraviolet initial scale A, where 
the flow is initiated. The latter can be identified with the 
ultraviolet momentum cutoff in the DSE framework. 

The work is organised as follows. In Section [IT] 
we discuss the FRG framework for (/(^-theory, includ¬ 
ing its properties and limits of the current simple 
approximation. In Section [llT| the formulation of the FRG 
at finite temperature and volume is introduced in a self- 
contained way, with a special emphasis on the fate of 
long range correlations and the (non-)existence of con¬ 
densates. In Section |IV| we contrast thermodynamical 
quantities in an infinite and finite volume, and their FRG 
evolution. Special attention is paid to the relation of the 
chosen cutoff procedure (regulator function) and the de¬ 
cay behaviour of thermal and finite volume effects for 
large cutoff scales. In Section |V| we extend the decay 
analysis to the full dynamics. We also present numer¬ 
ical results for the pressure, the dynamical mass and 
(/(^-coupling for finite volume and temperature adapted 
regulators. 


II. (^'‘-THEORY AND THE FUNCTIONAL 
RENORMALISATION GROUP 


We study finite volume effects at the simple example of 
a one-component scalar field theory. Its classical action 
is given by 


5 [</.] = f 






p= 
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with the standard kinetic term and a potential 

= ^{P - - ^{Po - : Po = max(0, k) . (2) 

The potential comprises the mass term and the (j)^- 
potential with coupling A. The minimum of the poten¬ 
tial is at po = max(0, k). The potential only depends on 
(j}^ , thus the theory is Z 2 -invariant under </> -(f). The 

parameterisation of the potential in Q has the normali¬ 
sation l/'^*(po) = 0 at the minimum pg. 


A. Functional renormalisation group 

The functional renormalisation group (FRG) is a 
functional continuum approach that includes non- 
perturbative effects. In this Section we highlight those 
aspects which are relevant for this work. 

The idea of the FRG is based on Wilson’s idea of inte¬ 
grating out fluctuations momentum-shell-wise. In prac¬ 
tice, this is done by starting at a microscopic scale A, 
where the theory is defined via its classical action S'[(^]. 

In this work, we consider a field theory for a (one- 
component) scalar (f>. Fluctuations from larger distance 
scales with respect to the microscopic scale, i.e. energies 
lower than the renormalisation group (RG) infrared scale 
k, are suppressed by adding a regulator term, Rk{p), 
which leaves the ultraviolet unchanged but serves as an 
effective infrared mass ~ . This allows to integrate fluc¬ 

tuations which are of the order of this mass. By variation 
of the scale from microscopic to macroscopic values fluc¬ 
tuations from all scales can be taken into account succes¬ 
sively. Within this procedure, one smoothly interpolates 
between the classical action at large cutoff scales k and 
the full quantum effective action r[(/)] at vanishing cutoff 
scale, k = 0. This is achieved in terms of an integro- 
differential equation, the Wetterich equation, [T5] . 

dtTk[(l)] = 2 + 

with the standard abbreviation = f d^g/(27r)'^, and 
t = Infc/A with some reference scale A. The schematic 
notation on the right hand side of <§ stands for inte¬ 
gration of the full, non-perturbative diagonal part of the 
propagator I/(r(^^ + Rk)(p,q) and the derivative of the 
regulator Rk{q) over the loop momentum q. The initial 
condition for (|^ is given by rA[<()] = 5” [(/)]. For gen¬ 
eral k, the effective action rfc[())] describes the full quan¬ 
tum effective action for momentum scales larger than fc, 
and lacks the quantum fluctuations of momentum scale 
smaller than k. In the limit fc ->• 0 it turns into the full 
quantum effective action F[^]. 

The regulator implements the shell-wise integration of 
fluctuations described above, i.e. suppressing the infrared 
while leaving the ultraviolet unmodified. Within these 
restrictions it can be chosen freely. In the following, we 



FIG. 1: Different regulator functions, see @, § and 0 - 


study various common choices like the exponential regu¬ 
lator, 

where the parameter m controls the sharpness of sup¬ 
pression above the scale k, the flat regulator [TB] , 

Rt\q) = {k^-q")Q{k^-q^) , ( 5 ) 

which is optimised within the lowest order of a derivative 
expansion, see pRtiTH] . and the sharp regulator, 

The sharp regulator implements the standard ultraviolet 
momentum regularisation within Dyson-Schwinger equa¬ 
tions, as well as in perturbation theory. There, loop mo¬ 
menta q^ are cut at an ultraviolet scale A with q'^ < A^. 
Integrating the flow with the sharp regulator from an 
ultraviolet initial scale k = A to k = 0 precisely leads 
to such loop integrals: in this case the flow gives the 
standard BPHZ-type renormalisation. Other regulators 
give a generalised BPHZ-type renormalisation, see |18j . 
in particular the regulators Q do not spoil the analytic 
properties of the integrands for finite momenta. As we 
shall see, the non-analyticity of the flat regulator is al¬ 
ready causing large intricacies for the flow of thermody¬ 
namical observables. These intricacies get enhanced by 
the stronger non-analyticity of the sharp regulator. 

The regulators 0, 0, 0 are shown in Fig. 0 Note 
that the regulator does not need to depend on both, fre¬ 
quencies and spatial momenta. In fact, in the context 
of thermal field theory it is convenient to regulate spa¬ 
tial momenta only, q^ ->• q^ in the regulators defined in 
0-0, because in this case Matsubara frequency sums 
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can often be performed analytically. We will elaborate 
on this in section IV where we focus also on the suitabil¬ 
ity of regulators from additional constraints imposed by 
finite temperature or volume. 


B. Local potential approximation and flow of the 
effective potential 

For large cutoff scales k A, the effective action 
Ffc tends towards the classical action S in 0 with A- 
dependent parameters. In particular the effective poten¬ 
tial Vk{p) tends towards the (bare) classical potential, 
Va = Ki- Here, the effective potential 14 is defined as the 
effective action F^ evaluated for static (constant) fields 
Pc; 


Vk{Pc) 


^k[p=Pc] 

V4d 


with Vid 



(7) 


In Q the trivial volume factor is removed. The effective 
action also contains higher derivative terms which are 
dropped here in the spirit of a low energy derivative ex¬ 
pansion or local potential approximation (LPA). In this 
approximation, the full effective action reads 

Tk[<P] = f d^x(^^{d^cl)f + Vk{p)), ( 8 ) 

Within the LPA approximation (|^ , the fc-dependent full 
propagator is given by 


Gk{q,p) 


1 

+ ml{p) + Rk{q) ’ 


with the k- and field-dependent mass function 


Wfc(p) = Vkip) +‘2pVk (p) ■ 


(9) 


( 10 ) 


The primes indicate derivatives with respect to p. The 
effective potential Vk{p) can be expanded in powers of 
(p “ Po) about the equation of motion (EoM) for static 
fields. 

This leads to 

Vfc(p) =I4(po)- ^(po-k)^+ ^{p-nf 

+ I^(p-Po)", (12) 

with the relation between po and k from 0 and Va(po) = 
0. The first line comprises the normalisation of the 
effective potential and the fc-dependent counterpart of 
the classical potential including fluctuation effects for k 
and A, while the second line comprises the fluctuation- 


induced higher order scatterings. In this expansion the k- 
dependent analogues of the power-counting relevant clas¬ 
sical parameters A, n are determined by 

^k(^P0,k ~ ^/c) ~ ^k^PO.k^ ^ Afc = 14 (po,k) 5 (1^) 

where we made the ^-dependence explicit. The mass 
function ml(p), evaluated on the EoM, reads 


^k ~ ^A:(P0,fc) “ A/i;(3po,fc kik') 


(14) 


by using the definition of k in ([^. Vanishing k signals 
the phase transition. For large cutoff scales fc A, the 
effective potential tends towards the classical potential 
Vc\ given in 0. 

The flow equation for the effective potential is easily 
derived from 0 and reads 


Hfe(p) = ^ jG{q,p)Rk{q ), 


(15) 


with G{q,p) defined in 0 , and the dot indicates the di¬ 
mensionless derivative with respect to the RG scale, i.e. 
Vk = dtVk- The flow of the effective potential, (15) de¬ 


pends on V' and V via the mass function ml{p), (101 in 


the propagator. By taking the first and second derivative 
of (15) with respect to p in ([^ we find 


dpVk{p)=^ [ ^Gk{q,p)Rk{q) 

Jq 


2 Jq dp 


and 


= -^{ml{p))'J^Gl{q,p)Rk{q), (16) 

9p^fc(p) = [(wfc(p))'] J^Gl{q,p)dtRk{q) 

-^{mlip))" J Gl{q,p)Rk{q). (17) 


Eq. (16) and ( [T7| ) depend on first and second p-derivatives 
of the mass function, whose flows can be obtained by 
further derivatives of the flow equation 0 with respect 
to p. 


C. (;A^-approximation 


The partial differential equation (15) can be solved 
within various methods. Here, we resort to a Taylor ex¬ 
pansion about the minimum pg. Then, evaluated at po, 
([^-(@ are the first three of an infinite hierarchy of 
ordinary coupled differential equations for the couplings 
K, A, A „>3 defined in (12). At the initial scale A we have 
An,A = 0 for n > 3. Using the parameterisation (12) we get 
for the left hand side of the flow of the effective potential 
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([15]), 

dtV(p) =V(po) - po [V'{p) - V'ipo)] + \{po -k){p- po) 
+ ^{p- + E ^{p-PoT ^ 

2 2 ^3 n! 

(18) 

leading to 

V'{po) = -Xk, V(po) = X - P 0 X 3 . (19) 


with the solution for and Afc. 

In the current ^^-approximation the mass function and 
its derivatives are given by 

'm\{p) = Xki'ip-Kk) , ml' = 3Xk, ml'' = 0. (23) 

In summary, this leads to the closed system of flow equa¬ 
tions for fk, Kk and Afc, with 

fk = ^ J^Rk{q)Gk{q), 


The A 3 term in (191 signals the ^-dependence of the ex¬ 


pansion point and feeds into the flow of A. This exem¬ 
plifies the destabilising backreaction of the higher cou¬ 
plings on the flow of the lower couplings in the presence 
of a flowing expansion point, and can be avoided with 
a fc-independent expansion point, m- In the present 
work we drop the higher couplings also for fc < A, that is 
A „>3 s 0 and approximate 


1 


Vkip) ~ ’{po) + {p-Po) Vfc'(po) + - (p-po) Vl'ipo ), 

( 20 ) 


with (13). This is simply the first line in (12). The rapid 
convergence of this expansion in 0(A^)-models at inhnite 
volume in general dimensions has been well-studied for 
example within the computation of critical exponents, 
see e.g. PDII23) . This approximation also works well 
for low-energy effective models for QCD with mesonic 
degrees of freedom at vanishing density and finite tem¬ 
perature |19] . In turn, at finite density and finite tem¬ 
perature higher powers in p, further couplings and mo¬ 
mentum dependences are required for quantitative state¬ 
ments, for details see Ref. m- This concerns in par¬ 
ticular the existence and location of a critical end point 
(CEP) in the phase diagram. Then one either resorts to 
a ^-independent expansion point or relies on global tech¬ 
niques for solving partial differential equations, see e.g. 


Here, we are interested in structural results and we 
only keep the relevant parameters k and A for the sake 
of simplicity. Then the set of flow equations ( [T^ - 0 
is closed. In this approximation, we have acce ss t o the 
flowing free energy density fk, the mass m^, see ( [T^ , and 
the coupling Afc with 


fk = Vfe(Po) 


rfc[po] 

V4a 


( 21 ) 


kk = lf^Rk{q)Gl{q), 

Xk = 9XlJ^Rk{q)Gl{q). (24) 

Note that all flows are positive and therefore leading to 
decreasing functions fk,Kk and Afc. In summary this re¬ 
sults in an increasing mass. This reflects the property 
that bosonic flows are symmetry-restoring. Note also 
that the reduced system of equations for the mass and 
the coupling is already closed. The flow of ft only de¬ 
pends on TOfc = Afc|Kfc|. The above flows integrations over 
the four-momentum q: for finite extent in either tempo¬ 
ral or spatial direction the integration (partly) turns into 
summation over discrete modes. We elaborate on these 
modifications in section IIIII 


III. FINITE TEMPERATURE & VOLUME 


At finite temperature and/or in a finite volume, some 
of the space-time directions only have finite extent, 

Xo6[0,1/r], and/or Xie[0,Li]. (25) 

Note that in the current Euclidean formulation there is 
no technical difference between the spatial and temporal 
directions. One might as well interpret a finite extent in 
the temporal direction as a theory at vanishing tempera¬ 
ture and one compact spatial direction with Lq = 1/T. A 
finite extent in a given direction with periodic boundary 
conditions (for bosons, (j){x + L^) = (t>{x))) only allows for 
plane waves that are periodic under shifts x^ x^ + L^, 
that are exp{i2TTnx^jL^) with n € h. This entails that 
frequencies and momenta are discrete, 


and 


ml = Vl{po)+2poVl'{po), Xk^Vl'ipo). ( 22 ) 


The two parameters in ( |22[ ) can be determined from the 
coupled set of flow equations for dpVk{po) and dpVk(po). 
The flow of the free energy density fk in (21) only de¬ 


pends on Kk and Afc, and can be integrated separately 


Pp = 27rn— , with n eh and T = — . (26) 

Lp Lq 

In the case of finite temperature these are the Matsubara 
frequencies. Accordingly, spatial momentum integrations 
turn into sums over spatial Matsubara modes. The com- 
pactification of a coordinate x related to the momentum 
p to the interval a; € [ 0 , L] with periodic boundary con- 
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ditions leads to 



This is trivially extended to higher dimensions. 

In the case of finite temperature the temporal direction 
xq is compactihed to the range xq e [0, 1/T] with periodic 
boundary conditions for bosons. For fermions one has 
anti-periodic boundary conditions. For infinite spatial 
extent, momentum integrations have the form 

[ I{q^) ^ [l{{27rTnof + f), (28) 

-'9 Tlo = -O0-'? 

with the spatial momentum integration f^ = f d^(7/(27r)^, 
and a sum over the Matsubara modes 27rTno. 

In the case of finite temperature and finite spatial ex¬ 
tent L we are left with sums over discrete modes in both 
temporal and spatial directions, respectively, viz. 

Jq L<oo 


T 

1 ? 


E ^ 


n^=-oo 

M=0,1,2,3 



2 9 2 

0 2 + no + no 


(29) 


with ui, 712 and label the discrete spatial modes in 
the cubic volume of edge length L. 


A. Range of finite temperature & volume effects 

Finite extent, ( |25[ ), seemingly induces a dimensional 
reduction of the theory in the presence of finite temper¬ 
ature and/or finite volume: for large temperature/small 
volume only the zero mode of the sum contributes, leav¬ 
ing only the integration/sum over the other momentum 
directions. This dimensional reduction is well-known at 
finite temperature. There it is commonly formulated that 
a quantum field theory in 3 -r 1 (or other dimensions) is 
dimensional reduced to its 3-dimensional counterpart for 
sufficiently large temperature. However, strictly speak¬ 
ing this only holds for momentum scales pjT « 1. There, 
however, it is valid for all temperatures. In turn, for 
p/T » 1 the dimensional reduction does not take place. 
This is reflected in the fact that the renormalisation can 
be chosen to be temperature-independent. 

The co-existence and interaction of these two mo¬ 
mentum regimes for all temperatures and/or finite 
lengths has the direct consequence that we see both, 
an exponential suppression of thermal and/or finite 
volume effects with the physical mass scales of the 
theory as well as the standard polynomial suppression of 
large momentum modes already present in perturbation 
theory: 

Exponential decay: In finite temperature held theory 


thermal correlations decay exponentially with the mass 
gap TOgap of the theory, which can be formalised as 


lim 




T=0 


QT=0 


2l(^)exp(-(IS2;). (30) 


for correlation functions 


= (</(a:i)-"(/>(a:„))T, 


(31) 


in the absence of further scales. If the latter ones 
are present they lead to further dependencies on their 
dimensionless ratios in A. The function A{y) is a 
rational function of the argument y = irigi^plT, and the 
exponential suppression with ejq){-mgnplT) can be 
readily computed from one loop thermal perturbation 
theory. Evidently, similar expressions hold for the hnite 
volume correlations, with the identihcation T ->• 1/Li 


in (30). For potentially vanishing correlation functions 


the denominator of (30) should be substituted with 
\^n\ + the sake of definiteness. For the 

correlation functions studied in the present work this is 
not necessary. 

Polynomial decay: 


Note that (30) does not im¬ 


ply an exponential suppression of thermal corrections 
in momentum space with momenta proportional to 
exp(-const.p/T). In momentum space, correlation func¬ 
tions show a subleading momentum behaviour with 
temperature-dependent coefficients that is only sup¬ 
pressed by powers of p^. For the sake of simplicity, we 
consider the symmetric point p^ = p^ for all i = l,...,n 
and obtain 


lim 


OL -O 


T=0 


0^=0 


ocBl- 


(I) 


0 , 


(32) 


with a rational function B^pjT) which decays polyno- 
mially for pjT ->• oo. 

For the two-point correlation functions, the propaga¬ 
tor G{p), this is known as the Tan-contact term in e.g. 
the context of ultracold atoms P5I127) . Within the FRG 
these terms have been discussed in lilHlIlH] for Yang- 
Mills theory and ultracold atoms, respectively. In [Si] 
a detailed analysis of the diagrammatic origin of this be¬ 
haviour was provided. This also entails that the coeffi¬ 
cients of this decay come from thermal fluctuations and 
their values hinge crucially on the correct implementa¬ 
tion of the exponential decay in (32). Hence, even though 


these terms have a polynomial decay, the correct compu¬ 
tation of the prefactor Ct relies on the exponential sup¬ 


pression (30) with the physical mass scales of the theory. 


For the propagator we have 


lim 

B^O 


{0,p)-G^=^{0,p) 


G^=o(0,p) 


oc Ct I Tzr 




(33) 


with a temperature-dependent prefactor Gt (in the pres¬ 
ence of other scales), the Tan-contact. Importantly, the 
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Tan-contact is related to e.g. the equation of state as 
well as the density. Hence, it is expected that only an 
approach that is able to deal with both the exponential 
suppression and the polynomial terms enables us to com¬ 
pute these quantities. 


B. Flow equation at finite temperature & volume 


The flow equation for the effective potential takes the 
form 

T ■ 1 

^k{p)=^ E Gk{q,p)R{q), with q^ = 2TTn^—. 

(34) 

Note that the effects of finite volume are qualitatively 
similar in other dimensions (c? - 1) + 1. Thus, for the 
purpose of this work, we restrict ourselves to (3 + 1) di¬ 
mensions. In the following, we discuss results obtained 
in the ^^-approximation for flows at both finite temper¬ 
ature and finite volume. These results are obtained by 
integrating the flows (containing loop integrations (29l), 
(24), over the RG scale k from A to 0. 


In the current FRG setup the dimensional reduction 
discussed in the last section |III A| is apparent from the 
flow. To begin with, at finite temperature and fc/T « 1 
the flow only receives contributions from the Matsubara 
zero mode (for regulators that decay sufficiently fast for 
frequencies q^jT « 1). Consequently, in this regime the 
flow reads 


T 


Rk{qo = 0,|g|) 


Vk«T{p) 2 fq q'^ + Rk{qo = 0, \q\) + mk{p) 


(35) 


The prefactor T can be absorbed with an appropriate 
rescaling of 14 > P, leading to dimensionally reduced flows 
for 


V,, = VIT, p,, = pIT. 


(36) 


The new variables ( |36[ ) have the momentum dimensions 
of the 3-dimensional potential and field-squared. More¬ 
over, (35) takes the form and the dimensional properties 


of a (3-l-0)-dimensional flow for the potential. 

With the same line of arguments it can be easily 
seen that a finite volume setup does not support a non¬ 
vanishing condensate, effectively reducing to quantum 
mechanics for a large correlation length: if the cutoff 
scale k is far smaller than the inverse size of the system, 
k « 1/L, only the spatial zero mode contributes to the 
flow, effectively reducing it to a 0 + 1-dimensional flow 
with only the frequency integral to be performed. This 
mimics a quantum mechanical flows, to wit 


Vk«iiL{p) ^ f 


dqo 


Rk{qo,\q\ = 0) 


27r q^ + Rk{qo, I?! = 0) + ml{p) ' 

(37) 


As in (35) the prefactor l/L'^ can be absorbed with a 


similar rescaling of 14 and p similarly to (36), with 


V,^ = VL\ p,, = pL\ 


(38) 


This renders (37) in a form that is identical to that of 


the quantum-mechanical (O-l-l)-dimensional flow for the 
potential. The latter flow does not admit a non-trivial 
minimum for the potential: there is no spontaneous sym¬ 
metry breaking in quantum mechanics. 

Finally, at finite temperature and for k « min(l/L,T), 
the system reduces to the zero mode. 


^«min(l/J^,T) (7*) 


T 


k^ + mf.{p) 


(39) 


where, without loss of generality, we have used that 
Rk{0) = k^ and Rk{0) = 2fc^. Again this can be mapped 
into the standard flow by a rescaling of I4 and p similarly 


to 


138) with 


Vo, = VL^lT, p,, = pL^lT, (40) 


i.e. effectively taking T/L^ 1. This is the flow of a 
(O-l-O)-dimensional field theory, the generating functional 
being the one dimensional integral 

Zod,k{J) = f {-Hd,A(p) - + J 4>od} , (41) 

where A « min(l/L, T). This trivial flow has been stud¬ 
ied in e.g. Ref. [29] : it does not admit a non-trivial min¬ 
imum for the respective effective action Fq^[</>]. 


C. Condensates 


In a finite volume no long range spatial fluctuations 
are permitted. However, the latter ones are required for 
phases with spontaneous symmetry breaking and non¬ 
vanishing condensates. Only in one or more spatial direc¬ 
tions and with infinite extent (at vanishing temperature) 
non-vanishing condensates survive for the present Ising 
universality class model (O(I)-symmetry). 

The disappearance of condensates in a finite volume is 
easily seen in 0(Af)-models for > I due to the pres¬ 
ence of massless Goldstone modes with mass functions 
eip) = ^kip)’ where 6 are the A^-1 Goldstone modes. 
Here, emphasising the differences to the Ising-type case 
N = 1, we briefly discuss this case. 

The flow of the minimum can be derived from the t- 
derivative of the equation of motion, VKpq) = 0 for po > 
0. Its flow reads 


^iVkiPo) - V'{po) + pqV''{pq) - 0 , 


(42) 


where the flow I4' is given in (16). The masses of the 
Goldstone modes vanish in the broken phase, i.e. at a 


non-trivial minimum po > 0- By resolving (42) for po we 
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arrive at 


Po = - 


V'{p^) 

V"ipo) 


■llMq) [(3 + Po^)Gl^{q) + iN-l)Glg{q) 


1 T \ ^ + Pot 


fc2 L3 


+ (iV-l) 


(43) 


Still, the vanishing of the condensate at finite volume 
and temperature is easily seen. For the sake of complete¬ 
ness we briefly present the related argument. Assume for 
the moment that w|^q(po) > 0- Then, for k > 0 the po¬ 
tential has a non-convex regime with to|(p < p^,) < 0 and 
the turning point pw with m^{pw) = 0. For this regime 
the relative strength of the free energy AI 4 = 14 (p = 
0 ) - Fk is governed by 

^yk = ^ J^Rkiq) {Gk{p,q)-Gk{q)) , (45) 


The first equality in ( |43[ ) provides the flow of po in the 
local potential approximation in an 0(A)-model. The 
first term in the second line is the radial contribution 
while the second term comes from the Goldstone modes 
with the massless propagator Gg. In d < 2 (infinite range) 
dimensions and N > 1 the latter term does not tend to 
zero in the limit A: -> 0 while the first term does with 
powers of The Goldstone term is positive, and 

hence the condensate is driven to zero with k ->■ 0, that 
is t ^ - 00 . This reflects the Mermin-Wagner theorem. 


In a finite volume and at finite temperature the flow 
reduces to the third line for k « min(T, 1/L). Again, the 
singularity in the flow is apparent. As a consequence, po 
is driven to zero with fc -> 0. This is also seen in Yukawa- 
type models such as the linear quark-meson models and 
atom-condensate models for ultracold atoms, see e.g. P- 

EJ. 


Note that this argument does not apply to N = 1 with 
the discrete Z 2 -symmetry. The 0(1) model (Ising uni¬ 
versality class) has no Goldstone mode and the vanishing 
of the condensate for d < 2 dimensions cannot be read-off 
directly from the flow of the condensate, (43). In partic¬ 
ular, the flow in (43) reduces at finite temperature and 
volume to 


Po 


. 2 T 3 + po^ 

L3 (fc2 + 2poA2)2 ’ 


(44) 


for k « min(r, 1/L). This relates to the fact that the 
Ising universality class, N = 1, shows a phase transition 
in d = 2. Below two dimensions, d < 2, it does not show 
a phase transition. 

We remark that this structure is not seen in a Taylor 
expansion about the flowing minimum po- For example, 
in the (/)^-approximation the condensate does not neces¬ 
sarily vanish at vanishing cutoff scale: for a finite mass 
mf. the flow of all couplings, including po and the free 
energy F^, can vanish with powers of k^jirif.. Indeed, the 
non-vanishing minimum stays at higher orders at this 
expansion too. This hints at a failure of the polynomial 
expansion about pq. Note, however, that an expansion 
about p = 0 works in dimensions d < 2, and converges 
rather rapidly with the full solution of the flow equation 
for the potential. This suggests an expansion about the 
flowing minimum po for fc > 1/L and one about p = 0 for 
k « 1/L. 


where p < p^, is a, /c-independent field value in the non- 
convex regime with m^(p < p^,) < 0, potentially leading 
to large flows close to the singularity triggered by k^ + 
'rnl{p) = 0. For the current purpose it is sufficient to use 
the estimate 


1 


< 


p^ + Rk{p^) ' + Rk{p^) + ml{p < p^) 

and, hence, we estimate 

AI 4 > ^ ^ ^ 

Z Jq \ 

1 r R .(„\ I 

In the limit fc/[min(l/L,T)] ->• 0 we arrive at 

T 



AVk 


(46) 


L3 ■ 


(47) 


(48) 


Evidently, the flow of AI 4 does not vanish for fc ^ 0 
and hence the difference turns negative: the minimum is 
at p = 0. Note that an iteration of this argument also 
excludes the existence of several minima and, therefore, 
guarantees a smooth vanishing of the minimum po with 
fc^O. 

For vanishing temperature or more generally infinite 
range dimensions d < 2 the above argument can be gen¬ 
eralised by discussion of the flow of the turning point 
Pw of the radial mass function rn\{p) with m\{pw) = 0. 
We have p^ < po for the 0(l)-model and p^, = po in the 
0{N > l)-model. The full analysis is beyond the scope 
of the present work and hence discussed elsewhere. 


IV. THERMODYNAMICS IN A FINITE 
VOLUME 


One of the main goals of this work is the study of finite 
volume effects in the pressure, 

p = -dFjdV , (49) 

with the free energy F and the spatial volume V = L^. 
The free energy F relates to the effective action r[(/)EoM] 

















k 



FIG. 2: Spatial volume dependence of the flow of the free 
energy density /* for a coupling Aa = 0.5 and an exponential 
regulator with m = 1, cf. eq. 


on the solution of the EoM, ^eom- It can be normalised 
(shifted) such that it vanishes in the combined infinite 
volume and zero temperature limit, that is 




(50) 


Eq. (50) is the free energy related to thermal and finite 


volume fluctuations, and the pressure (49) with (50) is 


the combined finite volume and thermal pressure. In the 
following we use the free energy density fk already in¬ 
troduced in (21) in the cutoff scale-dependent case. We 
have 


„ _ ^ _ rfc[(/)Eo: 

■ V ^ 


with /3 = - , (51) 


by extracting the trivial spatial volume factor in the 
static free energy. Eq. (50) and (51) entail that the pres¬ 


sure pk comprises the fluctuations induced from finite 
temperature and volume. The difference of free energies 
densities fk at finite and infinite volume relates to the 
Casimir force, for a ERG computation see Ref. [30] . 

The flow of the pressure is then given by the one of 
the free energy density and its volume or, equivalently, 
length derivative. 


Pfc(r) = -/fc(r,L)-| ^^^^^’^\ (52) 

where we have used that V = L^. The flow of the free 
energy density fk is given by 

h{T. L) = \Y Gl'\q)Rk{q) Gl’-^{q)Rk{q) , 

(53) 


where the sign Y,q stands for a sum over frequencies and 
discrete spatial momenta. In Eig. [^the flow of the free 
energy density, (53), is shown for different volumina. 

In (53) the first term is the finite temperature flow 



FIG. 3: Flow of the Stefan-Boltzmann pressure psB,k = 
-fsB.kiT, oo) for different regulator, and the full pressure for 
in an infinite volume, with A a = 0.5 for the full result. 


involving the thermal and/or finite volume propagator 
G^'^{q), whereas the second term is the normalisation in 
terms of the vacuum propagator (^q'^ Moreover, 

in the present LPA approximation the flow of the free 
energy density fk at a given cutoff scale k is simply given 
by the difference of the flows of the effective potential, 

h = (54) 

cf. section [ll] The flow equation for the infinite volume 
pressure is given by that of the free energy density, 

Pfc(T) = - hm /fc(T,L), (55) 

L^oo 


as the L-derivative of fk vanishes faster than 1/L. The 
flow of the thermal pressure pk, (55), is shown for differ¬ 
ent regulators in Fig. [^ 

The influence of the second term in a finite volume 
can readily be discussed qualitatively. It is already well- 
known from purely thermal flows that /^ > 0 for k < 
min(r, 1/L). It is this cutoff regime where the main con¬ 
tribution to the pressure is collected during the flow, lead¬ 
ing to a positive pressure, see also Fig.[^ For larger cutoff 
scales the flow switched sign, /^ < 0 for k < min(T, 1/L). 
As the volume and temperature pressure grows with de¬ 
creasing volume and/or increasing temperature the above 
properties imply 

^^<0 for fc < min(T, 1/L) , 
dL y > I 


>0 for k> min(T, 1/L). (56) 

oL 


In summary the integrated contributions from 
-1/3 Ldkfk decrease the pressure in a finite vol¬ 
ume, and the free energy density fk=o gives an upper 
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estimate for the pressure. 


A. Stefan—Boltzmann pressure for 
three-dimensional flat regnlator 


A simple but instructive example for the properties 
and definitions discussed above is the tree-level free en¬ 
ergy density fsB,k and pressure PsB,k, be. the Stefan- 
Boltzmann contributions, for two reasons. Firstly, it il¬ 
lustrates on an analytic level the computations and re¬ 
sults presented in the next Section|V] Secondly, and more 
importantly, it demonstrates constraints on the regula¬ 
tors emerging from finite temperature or volume in a 
clear way. 

First, we compute psB in the infinite volume limit with 
the three-dimensional flat regulator, see ([^. This regula¬ 
tor is often used as analytic computations can be pushed 
far. The free propagator in the presence of is 

simply given by 


- 9 ") + - fc') ■ (57) 


Qo + 


95 + 9 


With the propagator in (57) and = 2k^9{k‘^-q^), the 
flow for the infinite volume Stefan-Boltzmann pressure 
(55) can be given analytically, see e.g. [HI 131], 


Ps.AT) = - I ^ 


Jo 27r2 \ Jqo k"^ + ql) 


127r2 


(coftT-l) 


For large RG scales compared to the temperature this 
turns into 


kjT- 



(58) 


where qs = and ^. The flow ( [s^ and 

that for other regulators is shown in Fig. Eq. ( |58| ) en¬ 
tails that for large cutoff scales k/T 00 the flow decays 
exponentially with exp(-fc/T) in line with ( |30[ ). This 
is expected as the cutoff used in ( [^ does not affect the 
frequency sum, and reflects the fast decay of thermal fluc¬ 
tuations in the presence of large mass scales. Hence, the 
pressure at the initial scale A » T can be safely put to 
zero, Psb,a = 0. The integration over the RG scale k from 
00 to 0 gives the correct pressure pse of a gas of free scalar 
particles. 


Psi 


iT)= rJ^PsB,k{T) = 

J OS k, 


2rp4: 




(59) 


The Stefan-Boltzmann expression (58) is easily upgraded 


to the full flow. In the latter case we simply have to 
substitute k^ ^ k^ + rnf. in the denominators in the first 


line of (58) and arrive at 


Pk{T)=- 


127r2 


1 + 


T2 


coth 


2T 


\ 


1 + 


mt 


127r2 




0 2 

fe2 / 


(60) 


The term in the first line of (60) is decaying exponen¬ 
tially for large cutoff scales fc/T ^ 00 also in the presence 
of non-vanishing mass function m^. In turn, the term 

in the second line is proportional to (m^) - : 

the exponential decay of the full flow in ( [60| ) hinges on 
the exponential decay of the difference of the mass func¬ 
tions at finite temperature and vanishing temperature. 
The flow in (60) already encodes the one of differences of 
the effective potential with mk ->• mk{p)- Evidently, the 

exponential decay for {rr{^{p)') - (m^=°(p)) is encoded 
in the flow. This simple but instructive example already 
illustrates the main features of the thermal and finite vol¬ 
ume properties under investigation: the flow reflects the 
thermal decay with large mass scales in case that this is 
also present for the flowing couplings. 

We now repeat this computation with the three- 
dimensional flat regulator in the general case with a fi¬ 
nite volume and at finite temperature. For the sake of 
simplicity we first drop the 9i,/fc-term in (52). The three- 
dimensional spatial integral turns into sums, and we ar¬ 
rive at the Stefan-Boltzmann expression for the flow of 
the free energy density. 


2L3 Zo 


q2<fe2 


2T 


khl 

k'^ + ql 



in a splitting similar to that in (60). In (61) we have 



with h = (ni,77,2, ns). 


(62) 


Results for different regulators are shown in Fig. 

The first term in parenthesis in the second line of (@ 
decays exponentially with exp(-fc/T) for large kjT, but 
does not decay exponentially with k L for fixed k/T. This 
signals the failure of the three-dimensional flat regulator 
to reflect the exponential decay of finite volume fluctu¬ 
ations in the presence of large mass scales in a similar 
decay with the cutoff scale. 

For a direct access to the finite volume decay we take 
the limit T ->• 0. Then, the first term in parenthesis 
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FIG. 4: Flow of the Stefan-Boltzmann free energy density 
-fsB,k for different regulators, and and the full free energy 
density, -fk, for in an finite volume with relative lengths 
LT = 1 and Aa = 0.5 for the full result. The colour coding 
and normalisation is the same as in Fig. 


vanishes and we are left with the second term. This limit 
gives the flow of the free energy density in a finite volume 
normalised with that in an infinite volume. It is expected 
to decay exponentially with fc L -> oo similar to a thermal 
flow. Instead, the combined flow fluctuates about zero 
with a rising enveloping amplitude proportional to . 

The sum in the second term simply counts the number 
of momentum points in a three-dimensional sphere with 
radius k. Within the numerical precision of the present 
work we find for a d-dimensional sphere, 

.jr'i/2 / kl \ 

£rF(H)(2^) 


where = (27r)^/L^(nf + ••• + n^) is the d-dimensional 
discrete momentum squared. For d = 1 we have trivially 
di = 0. Note that while the sub-leading behaviour in (63) 


has been proven for d > 4, it is still subject to ongoing 
research for d = 2,3, with d = 2 being the Gaufi circle 
problem. The best estimates we found in the literature 
are 6»3 < 21/16 [ 33 ] and 1/2 < 02 < 131/208 [ 31 ]. In Ap¬ 
pendix for illustration, we have performed a similar 
computation for a regulator that achieves a flat propa¬ 
gator in a spatial momentum box with box lengths 2 k. 
Such a regulator allows for fully analytic computations, 
and is even more adapted to the symmetries of the spatial 
box. 


With (63) the Stefan-Boltzmann flow /sB,fe in a fi¬ 
nite volume diverges proportional to k'^ for the three- 
dimensional flat regulator. This asks for a non-vanishing 
initial condition /sb,a at the initial cutoff scale A: = A in 
order to guarantee that /sB,fc=o is the tree level free en¬ 
ergy density. Corrections of the thermodynamics due to 



FIG. 5: Flow of the pressure PsB.k for different regulators in a 
finite volume with relative lengths LT = 1. The colour coding 
and normalisation is the same as in Fig. 


the full dynamics are encoded in the fluctuating part 

/ \ 


fk fsB,k °'- 


127r2 


1 


1 






k 


E ■ 

q^<k'^ 


127r2 


1 + 




(64) 


where we have dropped the exponentially decaying parts. 
The term on the right hand side of the first line in (64) is 
proportional to k"^ The mass dif¬ 

ference tends to a constant, albeit small (up to logarith¬ 
mic corrections), for large scales. In turn, the second line 
is proportional to the difference of three-dimensional in¬ 
tegral and three-dimensional sum, that is the sub-leading 
term in (63). Hence, its total scaling is proportional to 


The above analysis also entails that for the present 
regulator, the second term in the flow of the 

pressnre, ( [5^ , a priori cannot be dropped even in the 
infinite volume limit. It can be easily computed from 
(61), and reads 


L fsB,k 


■ k'^ 1 A;^ ^ 2 

-/sB,fc + “ 3 ^ “ 9 ) , 

(65) 


where the second term is vanishing everywhere except for 
kf = (27r/L)^(n3 + + n|), simply collecting the contri¬ 

butions at the discrete momentum values. Using both 
terms, (61) and (65), in the flow equation for the pres- 
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sure, (52), we arrive at the simple expression 


PSB,fc — 


127r2 3 L3 


^-^coth^ E S{k^-f)- (66) 


neZ^ 


The flow of the Stefan-Boltzmann pressure for the three- 
dimensional flat regulator is given by a fc^-term and a 
i5-function contribution that peaks at the discrete spatial 
momenta. For the full pressure flow we can repeat the 
same analysis as for the flow of the free energy density 
leading to similar /c^-corrections. 


In Fig. we show the flow of the free energy density 
for the three-dimensional flat regulator in comparison to 
their Stefan-Boltzmann counterparts. Again, the differ¬ 
ences are negligible for large cutoff scales as rn'llk'^ « 1. 
This limit is usually achieved in theories or models with 
small ultraviolet fluctuations in the validity range of the 
theories. Note that in the present ())‘^-theory the Landau 
pole is hit at a large cutoff ^Landau, and one has to safely 
stay below this scale. We also emphasise that both, the 
flow of the free energy density, -fk in (61), and that of 
the pressure, pk in (66), tend towards (58) in the infinite 
volume limit. For T ->■ 0 and L ->■ oo both flows vanish. 

The simple example of thermodynamical observables 
studied in the present Section has already taught us an 
important lesson: flows with non-analytic regulators do 
not reflect the physical exponential decay with large mass 
scales in the decay with the cutoff scale. While the latter 
property is not necessary in order to guarantee the former 
physical one, it is potentially cumbersome to achieve in 
approximations. 


B. Stefan-Boltzmann pressure for four-dimensional 
flat regulator 

For the four-dimensional flat regulator, , the tree- 
level propagator reads 

- 9^) + ^^(9' - fc"). (67) 

with the momentum squared only taking discrete val¬ 
ues 




rineZ, 


( 68 ) 


and n = ( 711 , 712 , 77 . 3 ). We have chosen all lengths to be 
equal, Li = L. With the propagator (67) and 7?fc‘“(g) = 
2k^9{k^ - q^), the flow of the free energy density is given 
by 


fsB,k - .3 E f 




T 


E - 


k^ 




• (69) 


Similarity to the computation above for the three- 
dimensional flat regulator, we have to add the second 



FIG. 6: The difference of the flow of the free energy density 
contribution for the four-dimensional flat regnlator 


term in (52). This leads to the simple flow 


PsB,k — 


2 T 

327r2 “ 3 Z 3 




(70) 


see Fig. The spatial and four- dim ensio nal m omenta 
squared ^ and are defined in (62) and (68), respec¬ 
tively. 

Similar to the case of the three-dimensional flat regu¬ 
lator we conclude that the flow of the free action density 
is sensitive to finite volume effects for all cutoff scales for 
the four-dimensional flat regulator. With a similar anal¬ 
ysis as in the last Section |IV A| it follows that there are 
subleading terms proportional to ((ru^’^)^ - ( 771 ^’°°)^) 


and fc2(j7j7' Pqj. choice rn^lk^ « 1, being the ini- 


tial value in both finite and infinite volume as well as zero 
and non-zero temperature, at large cutoff scales the cor¬ 
rections are small. In Fig. we show the flow of the free 
energy density for the four-dimensional flat regulator in 
comparison to their Stefan-Boltzmann counterparts, as 
well as the Stefan-Boltzmann flows for other regulators. 
In Fig. we show the difference of the flow of the free 
energy density. This difference clearly shows the ultravi¬ 
olet sensitivity of the flow with respect to finite volume 
effects. 


In the infinite volume limit both, (69) and (70) tend 


towards the infinite volume flow. For T ->• 0 and L ->■ 00 
both flows vanish. We also conclude that for both, 
the three-dimensional and four-dimensional flat regula¬ 
tors, the thermodynamics is approximately accessible for 
m^/A2 « 1 at the expense of a trivial non-vanishing ini¬ 
tial condition /a = /sb,a and pA = Psb,a- However, since 
the flow of fk and pk only depends on the couplings, the 
Stefan-Boltzmann terms and the A/^, Apk can be inte¬ 
grated separately: the total free energy density/pressure 
is simply given by 


rO dk 

Pk=o ~ Psb{T,L) + / —Apk, (71) 

JA k 

with the trivial Stefan-Boltzmann pressure Psb(F, L) in 
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a finite volume and the integrated flow of 

The findings of the present and the last section carry 
over to the theory in the presence of three- and four¬ 
dimensional sharp regulators, Indeed, in these cases 
the ultraviolet growth is very similar. Below (§ it has 
been discussed that the theory in the presence of sharp 
regulators is directly related to respective ultraviolet mo¬ 
mentum cutoff regularisations of DSEs in a finite volume 
and at finite temperature, [IMl]- Consequently the re¬ 
spective correlation functions miss the exponential de¬ 
cays unless one applies a volume and/or temperature de¬ 
pendent renormalisation procedure. Note that only for 
the present thermodynamical observables, pressure and 
free energy density, this shows up as a growth with the 
ultraviolet cutoff scale A, in higher order correlation func¬ 
tions the missing exponential decay is hidden in polyno- 
mially decaying terms. 

We add that these findings apply to all situations 
with a sharp momentum cut-off for discrete momentum 
modes. This includes in particular Landau level sums 
in the presence of magnetic fields. There it is well- 
known that the critical temperature, and therefore the 
phenomenon of (inverse) magnetic catalysis shows an ex¬ 
ponential sensitivity with respect to the scales of the 
problem. Hence it is particularly sensitive to the presence 
or absence of the exponential decay discussed above. For 
FRG-applications to the (inverse) magnetic catalysis see 
[fl5H39j , for other approaches with ultraviolet momentum 
cutoffs including DSEs and perturbation theory, see the 
recent reviews gnun]. 

We close with the remark that it is not simply the miss¬ 
ing differentiability of the regulator that spoils the decay 
with the cutoff scale. It originates in the non-analytic 
behaviour of the regulator at . In Appendix [b] we 

test the smooth modifications of the flat regulator sug¬ 
gested in [32] that are C°° &t at the expense of 

an essential singularity oc exp{-const./(l - p^/fc^)} for 
spatial momentum approaching k from below, p^ ^ k^. 
They show the same missing exponential decay in the 
flow, although the smoothening of the non-analyticity 
leads to smaller fluctuations at large cutoff scales, see 
Fig. [^ in Appendix This is not surprising as it is 
well-known from the thermal case that the exponential 
decay reflects the pole structure of the loop integrals in 
the complex plane if the Matsubara sum is rewritten as 
the (original) contour integral. This suggests that ana¬ 
lytic versions of the flat regulator have all the necessary 
properties for optimising flows at finite volume and tem¬ 
perature in the sense of the optimisation criterion in |18j . 


C. Pressure for general regulators 


The flow of the pressure for general regulators is based 
on that of the free energy density and its length deriva¬ 
tive. More generally the length-derivative of the flow of 


the effective action is given by 






— T Gh Rh G, 


dr' 


( 2 ) 


(q) 


dL 


(72) 


where the last term hits the L-dependences of the ver¬ 
tices. Eq. (72) can be inserted in the definition of fk, 


(51), to obtain the L-derivative of the general flow of the 
free energy density. In the present LPA-approximation 
the second term in the flow of the pressure, (52), reduces 
to 


3 dL ■' 


1 

+ — 
6 


2q dg2 


dml 


jG}^ Rh , 

(73) 


see Fig. This leaves us with the task of computing 
the L-dependence of the couplings. This is either done 
by performing computations at L and L + e and taking 
discrete derivatives or by using flow equations for the L- 
derivatives in a fixed volume, see Appendix 

We close the present Section with a brief summary 
of the results. In the Figs. IllSIlj the findings of this 
Section are illustrated. In Fig. [^ the Stefan-Boltzmann 
pressure in an infinite volume is shown for different reg¬ 
ulators and the full pressure for the four-dimensional flat 
regulator. One sees both the exponential decay with k for 
analytic regulators (in the frequency po) and the power 
law rise for non-analytic regulators (in the frequency). 
All flows are apparently smooth in contradistinction to 
the finite volume flows depicted in Figs. EE This is 
related to the fact that the chosen regulators with a go- 
dependence are functions of Qq + g^. This entails that 
the jumps for non-analytic regulators, that are related to 
the discrete Matsubara frequencies are smoothened out 
by the spatial momentum integration. Regulators of the 
form Rk(qo,q) = i?o.fe(9o) + Rs,k(q^) with non-analytic 
regulators i?o,fe l^ad to jumps similarly to that shown in 
Fig.0 

In Fig. the Stefan-Boltzmann pressure in a finite 
volume is shown for different regulators. Again one sees 
both the exponential decay with k for analytic regula¬ 
tors and the polynomial rise for non-analytic regulators. 
Additionally, one has discontinuities as described in (611, 

(69) for the flow of the free energy density and in (651, 

(70) for the flow of the pressure. The polynomial rise 
of the free energy density and the pressure with k^ and 
k^ in a finite volume for the three-dimensional and four¬ 
dimensional flat regulator, respectively, hint at similar, 
but reduced problems for the flow of the mass (or 
Kfc) in both cases. The A:-power counting of these flows 
has an additional factor k~‘^ naively leading to rising 
volume-dependence with k and respectively. If the 
naive power counting holds, some care is required with 
the initial conditions for the mass, as then it varies with 
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FIG. 7: Normalised mass difference Am^ at finite tempera¬ 
ture T = 0.25 in an infinite volume as defined in (76l. 


FIG. 8: Normalised (^^-coupling difference AA at finite tem¬ 
perature T = 0.25 in an infinite volume as defined in (76l. 


the volume and the temperature. 

For the coupling constant we expect a Ijk suppres¬ 
sion of the flow for the three-dimensional flat regulator 
and a fc*^-dependence for the four-dimensional flat regu¬ 
lator. This translates into a slow convergence of the flow 
towards the infinite volume flow in the first case, hence 
asking for a very large initial cutoff scale. In the lat¬ 
ter case this leads to a logarithmic dependence with k of 
the initial condition on the volume and temperature, and 
requires some care. 


V. COUPLINGS IN A FINITE VOLUME 

We have already discussed in the introduction and in 
more details in Section m that a quantitative grip on 
the thermal, volume and cutoff scale dependence is piv¬ 
otal for getting hold of in particular thermodynamical 
quantities and the equation of state. Moreover, the ana¬ 
lytic results in the last Section |IV] indicate the problems 
of non-analytic regulator choices such as the flat regula¬ 
tor and the sharp cutoff. 


A. Thermal & finite volume decay 


In the following, we discuss the quantitative behaviour 
of the flows of couplings and thermodynamical observ¬ 
ables in dependence of the choice of regulators. Despite 
the choice of regulator both finite temperature and vol¬ 
ume are pure infrared effects, their effects are damped 
above energy scales which are typically of the order of 
the first frequency either in temporal or spatial direction. 
In other words, for scales much larger than the volume 
or temperature the theory turns into the infinite volume 
vacuum theory. The details of the decay are physically 
important, for the decay behaviour in mass scales and 
momentum scales see (301, (32) and (33) for the impor¬ 
tant example of the Tan-contact. 


In methods where fluctuations at all scales are inte¬ 
grated out at once this usually happens around the typ¬ 
ical frequency or spatial momentum scale given by the 
first frequency 2ttT of 2tt jL. This feature does not nec¬ 
essarily hold for the momentum cutoff scale k flow. As 


has been already seen in Section IV for thermodynami¬ 
cal observables, this cutoff scale depends strongly on the 
shape of the regulator. More precisely, it depends on the 
way the infrared modes are suppressed: the sharper this 
suppression is implemented, the further in k one has to 
go to lose sensitivity to thermal or finite volume effects. 
In particular for non-analytic regulators such as the flat 
regulator and the sharp cutoff this loss of sensitivity is 
delayed to infinity. 

For the regulator shapes given in Q, ([^ and § reg¬ 
ulating either four-dimensionally or in spatial directions 
only, the thermal flow of the tree-level pressure PsB,k is 
given in Fig. In an infinite volume the problem of 
a ’’mode-counting” cutoff can be circumvented by regu¬ 
lating the spatial momenta only. In contrast, in a finite 
volume this generic property rules out sharp and flat reg¬ 
ulators for practical computations. We note that this is 
not cured by the inclusion of running masses and cou¬ 
plings, as also seen in Fig. [^anticipating results below. 

We conclude that the exponential thermal decay in 

g does not simply hold with rUgap k, even though 
has been observed for the three-dimensional flat reg¬ 
ulator, see (58). In general, one has to allow for a pref¬ 


actor ct with 


exp 


(-^^) ^exp(-CT^) , CT 6 [0,1]. 


(74) 


For regulators that only depend on spatial momenta and 
have an (effective) infrared mass i?fc(0) = k^ the upper 
limit, Ct = 1, is saturated: the sum over the Matsubara 
modes is not cut-off, and the thermal behaviour is un¬ 
touched. With the dominant mass scale fc for fc ^ oo this 
leads to the standard thermal suppression with ct = 1. 
Note that strictly speaking the upper limit is not well 
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defined as it depends on the mass scale introduced with 
the regulator. Only if the latter scale is identical with k 
the prefactor ct satisfies (1^). 


In turn, for regulators that are non-analytic in fre¬ 
quency space such as the four-dimensional flat regulator, 
the exponential suppression is lost and the lower bound 
is saturated, ct = 0. 

In a finite volume the decay towards the vacuum limit 
is induced by both, thermal and volume fluctuations. 
Hence, it leads to a generalisation of the thermal decay 
(30), (74). For (53), we parametrise the decay in the 


limit 


lim 

klT,k L-*oo 


\AOl’^\ oc 


fc^]exp^-CT(FT)^j , 


(75) 


where the function A now encodes thermal and volume 
effects, and 


^ _ — 

^^k - ^0,oo 


(76) 


AO^’^ is the difference of an observable at finite 
volume and temperature and inhnite volume and vanish¬ 
ing temperature, normalised by its flow at T = 1/L = 0. In 
Fig.j^and Fig.j^the normalised difference AOk is shown 
for the mass and the coupling difference respectively for 
finite temperature in an infinite volume. 


The dehnition (76) also applies to flows or diagram¬ 
matic kernels of flows, the threshold functions. Flow 
equations of thermodynamical observables and vertices 
depend on the cutoff derivative of the regulator and pow¬ 
ers of propagators. The decay of general flows can hence 
be measured by that of 


O 


T,L 


fGkiqrilkiq). 

Jo 


(77) 


In Fig. we show the decay of general flows for different 
regulators, precisely, we show 


_qO c 
k,n fc,r 

J^d-2{n-l) 


(78) 


with the definition (77) for n = 1,2,3,4. We normalise 
with respect to the canonical dimension In 

the presence of an exponential decay, that is ct 0 for 
all show the same decay as they can be 

generated by n - 1 derivatives with respect to m 

The thermal decay at infinite volume and the volume 
decay at vanishing temperatures are given by the respec¬ 
tive limits. For LT oo, (75) gives the thermal decay 


I 


(74) with Ct = ct(oo), that is 


lim 

klT^r> 


|AO 


oc H ^ —. oo 


k 

.T’ 


)exp(-CT^) , 


(79) 



FIG. 9: Flows given by (781 for n = 1,2,3,4 in a finite volnme 
for different regulators. 


see Fig. Fi g. This is the analogue of (30) with the 
substitution (74). In turn, for LT ^0, (75) describes 
the volume decay at vanishing temperature, and we have 


lim 

LT^O 


ct{LT) 

LT 


CL, 


(80) 


and the equivalent relation to the purely thermal decay, 
(79) for the pure finite volume case. 


lim |AC>fe(0, L)| oc H (oo, fc L) exp (-Ci fcL) . (81) 

Lk-*oo 


The function ct{LT) interpolates between the asymp¬ 
totic thermal scaling and the asymptotic finite volume 
scaling, and depends on the shape of the regulator, 
ct{LT) = CT{Rk',LT). It can be extracted in the limit 
of A: ->• oo by 

ct(LT) = -y lim log{|AOfc(T,L)|} . (82) 

/c k^oo 


In Fig. we study the approach of the right hand side 
of (82) towards ct for n = 1,2,3,4 and an exponential 
regulator, ([^, with m = 2, in both the infinite and finite 
volume case, in the latter case we use L = 1/T. The local 
minima indicate the local maxima and minima in the 
positive and negative envelopes of the decay, as seen in 
Fig.0 and Fig. 1^ for n = 1. We find that the asymptotic 
limits for k/T and/or kj (1/L) are reached quickly above 
the characteristic scale induced by ct- 


In Fig. we study the interpolation of Ct between the 
asymptotic forms at T = 0 and L = oo for the exponential 
regulators with n= 1,2. For this figure we have extracted 
Ct with the value of the minimum at maximal fc, in e.g. 
Fig. the minimum between 87rT and IOttT. We see 
that the interpolation regime at about LT rs 1 is rather 
small and the asymptotic values are approached rapidly. 
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scales k larger the temperature and volume as parametrised 
by (821 and the an exponential regulator, (j^, with m = 2. 


Thick (thin, lower) curves are in inhnite (finite with L = 1/T) 
volume. 
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FIG. 11: Interpolation of ct for an exponential regulator, (j^, 
with m = 2 from the asymptotic purely thermal decay to the 
asymptotic purely finite volume decay. 


B. Thermal & volume dependence of the couplings 


and 


The results obtained in the last Sections exemplified 
the thermal and finite volume decays at the example of 
the thermodynamic observables, the free energy density 
and the pressure. Both observables have power counting 
dimension four, and indeed their flows scale with for 
non-analytic regulators. This leads to non-trivial initial 
conditions for these observables at the initial large cutoff 
scale k = A. We expect only power-counting reduced 
problems for the flow of the vertices, in particular for the 
mass rnf. and the coupling Xk- If this reduction leads to 
sufficiently strong suppression in powers of 1/fc even non- 
analytic regulators would allow for a direct access to both 
the thermodynamics as well as the flow of the couplings 
without resorting to non-trivial initial conditions at least 
beyond leading order. 

As for the thermodynamics it is illustrative to simply 
consider the length derivative of the flows. For the ef¬ 


fective potential this can be directly read off from (72), 

(P, 


L 


dVk 

dL 


2q^dq2 


dmi 


I ^k kik • 


(83) 


The second term in (83) can be viewed as an improve¬ 
ment term. It does not dominate the flow. For the 
present purpose of a power counting analysis we drop 
it and concentrate on the first term. The flow of the 
effective potential reads for the three-dimensional and 
four-dimensional flat regulators. 


dvr^^ip) 

— dr~ 




coth ^ 


4 


1 + 


((p) 


/c2 


4 


1 + 


^(p) 




-5{k^-f), 

( 84 ) 




ip) 


dL 


r3 Aj 


1 + 


((P) 


5{e-q^), (85) 


fc2 


respectively, where the superscripts indicate the chosen 
regulator. The question of the sufficient decay that is at 
the heart of T and L-independent initial conditions can 
already be answered at one loop. Thus, we integrate the 
flows in (85) at one loop from /c = 0 to A. This simply 


removed the ^-function leading to 


\P) 


\P)- 


T rSdflat 


T r4dflat 


{p) 


ip) 


-lA y 

'6L3 


coth Aw 


1 + 


Kp) 


A2 


4 


1 + 


((p) 


A2 


o g2<A2 1 + 


flA^ 


l(p) 


( 86 ) 


A2 


where we have also integrated over the box size from in¬ 
finite length to L. The expressions for t he one loop cou¬ 
plings m\ and Aa are derived from (86) by two and four 


derivatives with respect to ()>, respectively. This involves 
terms t hat are proportional to the first and second deriva¬ 
tive of (86) with respect to respectively. Hence, we 
have additional suppressions with 1/A^ and 1/A'*, respec¬ 
tively. In turn, the sums in (86) grow proportional to A^ 
and A* on the average. Sweeping over a single discrete 
frequency and spatial momentum value leads to a jump 
of the effective potentials and with a height 

proportional to A and A® respectively. For the couplings 
this entails a suppression of the single jumps with 1/A 
and 1/A^ for m\, and l/A^ and l/A"* for Aa. 

In the following, we restrict ourselves to the efficient 
regulator of the form (|^ with to = 1. For this choice 
the transition to the vacuum theory happens at cutoff 
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FIG. 12: Thermal and finite volume dependence of the mass 
for large UV coupling and T/A = 1/40, and the exponential 
regulator, Q with m = 1. 


scales close to the first frequency, and the cutoff scale 
can indeed be identified with a mass scale rugap of the 
order of the cutoff scale k. 

In principle, there are three different scales in the prob¬ 
lem at hand: vacuum mass, temperature and volume. In 
order to study the interplay of temperature and finite 
volume only, we consider the massless theory in the vac¬ 
uum. Note that this implies in general rn^ + 0 for cutoff 
scale k + Q. The vanishing mass at fc -> 0 requires a fine 
tuning at fc = A: for a given initial coupling Aa we vary 
the initial value of the mass mv such that mk^o = 0. A 
generalisation to finite masses is straightforward but not 
the main purpose of this work. 

In summary, from the value mk=o we generate the ini¬ 
tial condition wa and Aa for mass and coupling flows, 
(24) with (29), to be solved in a self-consistent way. The 
functional behaviour of the coupling and, consequently, 
of the mass is significantly different for large and small 
values of Aa. For large initial values of the coupling there 
are strong changes in mass and coupling, i.e. the running 
is very strong. For smaller values of the coupling the run¬ 
ning is slowed down. This results in the fact that the cou¬ 
plings at A: = 0 are much closer to each other than at the 
initial scale. This behaviour is of course well know and 
explained by the fact, that the coupling appears quadrat- 
ically and cubically, respectively, as a prefactor in the 
flows of the mass and coupling. 

The effects of thermal fluctuations have been studied 
extensively in e.g. [JS] SI] for the (/i"*^-theory, for Yang- 
Mills theory this has been done in [SHEj: As finite temper¬ 
ature is an infrared modification of the theory, the theory 
is unchanged at scales larger than the typical tempera¬ 
ture scale, which is set by the first Matsubara frequency 
2ttT. Hence, above this scale the flow follows the vacuum 
flow. Below this scale, temperature induces additional 
fluctuations which accelerate the flow. For the mass this 
results in the generation of a thermal mass in addition to 
the vacuum mass. However, once this generated scale is 



FIG. 13: Thermal and finite volume dependence for the cou¬ 
pling for large UV coupling. Regulator, parameters and 
colour-coding are the same as in Fig. |12| 


reached the mass suppresses fluctuation in both flows of 
the mass and the coupling. As a consequence, the flows 
are frozen below this scale and the values of the coupling 
and mass saturate. 

In this work, we consider the thermal field theory in 
a finite volume. Similar to the compactification of the 
temporal direction due to non-zero temperature, the fi¬ 
nite edge length of the (cubic) 3d-box is a modification of 
the infrared. Therefore, its effect is qualitatively similar 
to thermal effects. 

Once the theory is probed at distances comparable to 
the volume fluctuations become stronger and generate an 
additional mass. At some point, however, the flows stop 
as fluctuations are suppressed by the mass. In general, 
it depends on the magnitude of the three scales L, T 
and uiQ , as already the largest scale effectively suppresses 
long range fluctuations and leads to an independence of 
further macroscopic details. Having tuned to mg = 0, the 
only relevant parameter that determines the behaviour 
of the theory is the ratio of spatial length L to temporal 
length 1/T, hence LT. Thus, as long as cutoff effects are 
correctly taken care of, i.e. A » 2ttT and A » 27r/L, see 
discussion above, different combinations of L and T give 
similar results as long as their product is kept constant. 
In other words, the effects of finite volume are always 
relative to the temperature scale and it is thus sufficient 
to study only one temperature for different lengths. 

In Fig. Fig. we show cutoff scale dependence 
of the mass and the coupling for large initial values of 
the coupling. In Fig. [T^and Fig. we show the cutoff 
scale dependence for small initial couplings. Note that 
the flows are damped by unlike the flow of the po¬ 
tential. Thus, the technical problem of ’’mode-counting” 
cutoff does affect the mass or coupling but becomes only 
relevant at the level of the potential. It is clearly visible 
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FIG. 14: Thermal and finite volume effects in the mass for 
small UV coupling. Regulator, parameters and colour-coding 
are the same as in Fig. |12|but the coupling from Fig. |15| 



FIG. 16: Temperature and volume dependence of the mass for 
a massless (lower curves) and massive (upper curves) vacuum 
theory with = 0-2 for Aa = 5 at the initial scale. 


that the effect of the finite volume is qualitatively equal 
to the effect of finite temperature for both the mass, cf. 
Fig. 


Fig. 


12 


15 


and Fig. and the coupling, cf. Fig. and 
Furthermore, the dependence on the finite vol¬ 


ume is monotonic in the sense that the larger the vol¬ 
ume the less additional flow is generated, i.e. the infinite 
volume limit is approached. Interestingly, the limit of 
infinite volume is approached very slowly, however, i.e. 
the value oi LT » 1. This is due to the fact that the 
finite volume affects three dimensions and is to be com¬ 
pared with the compactification of the temporal direction 
only. The mass flows for a massless and massive theory 
are shown in Fig. |16[ which exemplifies that the vacuum 
mass does not change the effect of finite temperature of 
volume qualitatively. This justifies our choice toq = 0. 



k/A 


FIG. 15: Thermal and finite volume effects in the coupling for 
small UV coupling. Regulator, parameters and colour-coding 
are the same as in Fig. 


In Fig.[^we show the dimensionless mass as a function 
of the value of the coupling at 2ttT. We choose this value 
to facilitate comparison with other approaches where the 
coupling is defined at that scale, e.g. perturbation theory. 
Again, finite volume effects strengthen fluctuations and 
the effective mass is therefore increased. Thus, the infi¬ 
nite volume limit for the mass is approached from above. 
This is similar to the well-known thermal behaviour, see 
Fig.[l8l 

Eventually, we consider the pressure. Note again that 
is can be evaluated trivially by integrating the flow equa¬ 
tion ([^ with the mass function computed above. In 
Fig. [^the flow of the pressure is shown. Integrating 
the flow for a trivial initial conditions gives the pres¬ 
sure. In Fig. the pressure is shown as a function 
of the value of the coupling at 2ttT and normalised to 
the Stefan-Boltzmann pressure. The infinite volume case 



FIG. 17: Masses as a function of the coupling at 2'kT for 
different volumes. 
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FIG. 18: Mass ratio of thermal and vacuum mass as a function 
of (fixed) vacuum mass over temperature. The ratio shows an 
exponential decay with exp-m(0)/r. 


FIG. 19: Pressure, p, of <()^-theory in an infinite volume (black 
curve) and a finite volume (curves below infinite volume re¬ 
sult) in comparison with the free energy, -fk=Q (curves above 
infinite volume result). 


agrees with results in the literature. We see that both, 
increasing the coupling and decreasing the volume leads 
to a smaller pressure. Moreover, in comparison to the 
free energy density the thermal pressure shows signifi¬ 
cantly less volume dependence. This reflects the fact that 
it accounts for thermal fluctuations and their volume- 
dependence. 


VI. CONCLUSIONS 

We have investigated finite volume effects for the ther¬ 
modynamics of a (j)^ theory within the functional renor¬ 
malisation group approach. We have focused on concep¬ 
tual properties of the regulator dependence of flows in 
finite volume compared to infinite volume. The applica¬ 
bility of regulators can be decreased or even spoiled by 
the non-analyticity of the regulator for momenta at the 
RG scale, = k^. We find that the thermal flow of the 
effective action exhibits oscillations along k. Their range 
becomes larger with the sharpness of the regulator. In the 
limit of a sharp or flat regulator these oscillations extend 
to infinity. For a massless theory the only scales present 
are the temperature and the spatial length. Hence, ef¬ 
fectively the theory only depends on one parameter, the 
ratio of lengths in spatial to temporal directions, LT, 
where L is the edge length of a cubic box, and 1/T is the 
length in temporal direction. 

The cutoff scale introduces an additional scale to the 
system and both length scales are visible as onset scales 
for thermal and finite volume effects respectively. For 
analytic regulators the onset scales are proportional to T 
and L, while the proportionality constant depends on the 
chosen regulator. Above the onset scales the flows tend 
exponentially towards the corresponding vacuum flows. 


In turn, for non-analytic regulators the exponential de¬ 
cay towards the vacuum flows is missing, which makes 
quantitative studies of thermal and finite volume effects 
cumbersome. 

As expected, both finite volume and thermal effects 
show similar dependences on the mass scales of the the¬ 
ory including the cutoff scale. This is directly seen in vol¬ 
ume dependence of the thermal mass already, where the 
mass is larger for smaller volumes. As a consequence, the 
infinite volume limit of the mass is reached from above. 

Furthermore, we have studied the pressure. The infi¬ 
nite volume limit is reached quicker for larger couplings, 
because correlations are transferred in a more efficient 
way and fluctuations are washed out before the length 
scale is reached. Therefore, for smaller couplings the 
effects of the finite length is stronger. This is again a 
manifestation that these are pure infrared effects. 

The structural aspects of the finite volume depen¬ 
dence studied in the present work are quite general. 
The present findings carry over straightforwardly to 
0{N) models with more than one field component, 
to phenomenological low energy systems such as the 
quark-meson model, and, to some extend to full QCD. 
In the latter case the gauge sector requires some special 
attention. This will be studied elsewhere. 
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FIG. 20: Flow of the Stefan-Boltzmann free energy fsB.k for 
the css regulators ( |B1[) in finite and infinite volume compared 
to the fiat regulator in the infinite volume. 


Appendix A: Pressure flows for a flat box cutoff 


For a fully analytical approach we choose the product 
cutoff 


= n , (ai) 


fi=0 


with 


^ M=0 y 


i-Yle{e-ql) 


fi =0 


(A2) 


The flow of / is then given by 


/3B,fc=- (2iVT + l) 


'Nl 


+ Nl + 


2A|y 


+ 3{2Nl + 1)(—+N^ + 


2iV| 


Nl 


+ Nt + 


2^j^ 


277r4 


(A3) 


with 


Nl = 


k L 


Mx = 


k 

2 TTf 


(A4) 


where the brackets indicate the largest integer number 
smal ler than the argument. The enveloping function of 
(A3) rising as k^. This exemplifies the problem related 


to non-analytic regulators. 


Appendix B: Compactly supported smooth regulator 


Here we test the decay for compactly supported 
smooth (css) regulators suggested in [32], 

2^css i 


RT{<1) = q^rriq^lk^), 


rriy) = 


expjcyg/ (l-yg) - 1 } 

exp{c?/V(l-?/'’)-l} 


e(l-y), (Bl) 


with the normalisation such that r“®(?/o) = 1- We 
choose the parameters b = 1 , c = 1 and yo = 1/2 leading 
to a standard exponential decay. In Fig. we show the 
flow of the Stefa n-Boltzmann free energy /sB,fe for the css 
regulators (Bl) in comparison to the flat regulator ([^. 
In both cases of infinite and finite volume /sB,fe does not 
show the exponential decay for large cutoff scales k. 


Appendix C: Flows for the L-derivatives of the 
couplings 


In the present ^"^-approximation to the LPA one has 
to compute m^(L) and Afe(L). More generally, in a poly¬ 
nomial expansion of the LPA, where 

A = (Ai, A 2 , A 3 ,...), with Ai =/t, A 2 = A, (Cl) 

there are additional couplings. Their flows are given 
schematically by 


1 

K = - 

2 Jq 


JGk{q,po) Rk{q) (C2) 


and 


K-kXn^i=l- f Gi^\q,Po)Rk{q), (C3) 

Z Jq 


where G^'^\q,p) = d^G{q,p), using (15). Taking the L- 
derivative (or V-derivative) of (C3), with 


a; = L5lA„ = 3V 


dXn 

dV 


(C4) 


leads to flow equations for the L-derivatives of A„, there¬ 
fore, 

dtN = -I f {2q^dg2-X',dx,)Gl{q,po)Rk{q) , (C 5 ) 
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and 


dt^'n - dtX'iXn+l + dtXiX'^^i 

= -lf^{2fdg.-y,dx,)G^-\q,po)Rk(q). (C6) 


Eq. (C71 requires sufficiently smooth regulators in the 
spatial momentum directions. Moreover, iterations 
of this procedure leads to flow equations of iVth L- 
derivatives of the couplings A in the basis of j = 0, N -1 
L-derivatives of A and their flows. 


Note tha t po = max(0, k) depends on Ai = k. 
Eqs. (C5),(C6) can be solved on the basis of a solution 
for Afe in a given volume V (a cube with box length L), 
and a regulator that guarantees the decay of volume fluc¬ 
tuations for A: ->• oo, that is 


lim A^ = 0 . (C7) 


The present derivations trivially extend to the situa¬ 
tion of different spatial lengths Li for i = 1,2,3. We add 
that the present procedure applies more generally to flows 
of derivatives with respect to external parameters such as 
volume, temperature, chemical potential, and couplings. 
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